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TECHNICAL PAPER 


GIBBS FREE ENERGY OF REACTIONS INVOLVING SiC, SI 3 N 4 , H 2 , AND H 2 0 
AS A FUNCTION OF TEMPERATURE AND PRESSURE 


INTRODUCTION 

Silicon carbide (SiC) and silicon nitride (Si 3 N 4 ) are being considered as materials for 
coatings 1 and structural components 2 in hot sections of advanced propulsion systems using hydrogen 
(H 2 ) as working fluid. Possible applications include coatings for turbine and nozzle components. 
Advanced propulsion systems, including nuclear thermal, will most likely employ a feed system 
composed of a turbine and pump. Hot hydrogen will turn the turbine which drives the pump. 3 

A multiplicity of propulsion system designs is under consideration, yet no single design has 
been explicitly defined as the engine which will carry future craft on interplanetary flight. 4 Near-term 
designs consider turbine inlet temperatures up to 521 °C 5 and pressures up to 130 bars. 6 Nozzle 
temperature may reach 2,427 °C in the chamber 7 and pressure may be as great as 92 bars. 6 7 Higher 
pressures and temperatures generally contribute to higher engine efficiency and better performance. 
Higher temperatures across the feed system allow higher overall engine temperatures, contributing 
to more efficient operation. Higher reactor chamber pressure enhances performance through greater 
thrust. 5 7 

Several investigators have demonstrated strength reduction and weight loss in SiC and Si 3 N 4 
exposed to flowing H 2 as a function of temperature. 1 8 “ 10 Variation in moisture content affects 
material degradation, and a high enough level of moisture in H 2 may retard degradation by the 
formation of a silica (Si0 2 ) scale. Still, at constant moisture level, chemical reactions are 
temperature dependent with more severe degradation generally occuring at higher temperatures. 9 10 

System total pressure affects reaction equilibria. Misra 2 theoretically examined the stability 
of various ceramic materials in hydrogen gas at different temperatures, total pressures, and P H2 o- 

The partial pressure of gaseous products can be divided into discrete regions depending on moisture 
content of the hydrogen gas. Region I corresponds to a log(P H2 /PH 2 o) greater than about 10 5 . In 

region II log(Pn 2 /PH 2 o) is between about 10 2 and 10 4 . Silicon carbide and silicon nitride are generally 
reactive in region I, and they are actively oxidized to silicon monoxide (SiO) in region II. 2 

For SiC at all temperatures and pressures, methane (CH 4 ) is the major gaseous product in 
region I. In region II, SiO predominates at higher temperatures and lower pressures, while CH 4 
predominates at lower temperatures and higher pressures. 2 

In reactions of Si 3 N 4 with H 2 , nitrogen gas (N 2 ) predominates at higher temperatures and 
lower pressures in both regions I and II. As temperature decreases and pressure increases, SiFLj is 
favored slightly over SiO as the gaseous product. 2 

A relatively wet atmosphere (designated region III 1 2 10 ) is required to form a protective Si0 2 
layer on either SiC or Si 3 N 4 . Misra 2 has calculated that the moisture level required increases as total 


system pressure increases. It decreases slightly for SiC and increases for Si 3 N 4 as temperature 
increases. 

Reactions relevant to the decomposition of SiC and Si 3 N 4 in H 2 and H 2 /H 2 O are therefore 


SiC+2H 2 = Si+CR* (1) 

SiC+H 2 0+H 2 = SiO+CHt (2) 

Si 3 N 4 = 3Si+2N 2 (3) 

Si 3 N 4 +6H 2 = 3SiH 4 +2N 2 (4) 

2SiC+2H 2 0+2H 2 = Si+Si0 2 +2CH4 (5) 

Si 3 N4+2H 2 0+3H2 = Si+Si02+SiH4+2NH 3 +N 2 (6) 


Table 1 lists these reactions according to the conditions under which they are most likely to occur. 

Since service conditions may vary in temperature and pressure for SiC and Si 3 N 4 as structural 
materials and coatings, and degradation modes may vary with temperature and pressure, a fuller 
understanding of reactions operating at various temperatures and pressures is necessary. Flowing 
hydrogen gas creates a situation in which equilibrium is never reached. Fresh reactant is introduced 
continuously into the chemical reaction. Examination of the conditions under which SiC and S^N 4 
might preferentially decompose in H 2 is instructive. 


THEORY 


Reactions proceed spontaneously where differences between Gibbs free energy of products 
and reactants is negative. If AG for a reaction is positive, the reverse of the reaction proceeds. 
Change in Gibbs free energy is a function of temperature and pressure such that 

dG = ( 5G/5p) T dp+(Sg/$T) p dT . (7) 

Since (5 G/8p)r = V and ( 5G/8T) P = -S, equation (7) becomes 

dG = Vdp-SdT . (8) 


For molar quantities 


dGm = Vm dp-Sm dT , (9) 

where the subscript m denotes molar quantity. The integral taken over the ranges of temperature and 
pressure of interest is the total change in Gibbs free energy. 
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The change in molar volume of solids and liquids is negligible over a broad range of pressures. 
For practical purposes it can be considered constant. The molar volume of gases changes 
significantly with pressure. For an ideal gas 


Vm = RT/p . (10) 

Computation of molar entropy is similar for solids, liquids, and gases. It is 

Sm = Cpm/T (11) 

for solids and liquids, and 

Sm = CpmJT-R/p (12) 

for ideal gases. 11 Cpm is the molar heat capacity of the substance at constant pressure. 

The total derivative of molar Gibbs free energy with respect to temperature and pressure is 

then 


dGm(T,p) = Vm dp-Cpm/T dT 


for solids and liquids, and 


dGm(T,p ) = RT/p dp-Cpm/T dT+R/p 
= R(T+l)/p dp-Cpm/T dT 


(13) 


(14) 


for ideal gases. 

The Gibbs free energy of a substance at any temperature or pressure is found by integrating 
over the ranges of temperature and pressure of interest. If Gibbs free energy is determined with 
reference to the standard state (T = 298.15 K, p = 1 bar), it is the change in Gibbs free energy (AG) 
that is calculated. Integration gives 

AG(T,p) = AG°+Vm(p-pO)-Cpm ln(TIT0) (15) 

for solids and liquids, and 

A G(T,p) = AG°+R(T+l) ln(p/pO)-Cpm ln(T/T0) (16) 

for ideal gases. Where A G° is the change in Gibbs free energy of the substance in the standard 
state. Heat capacity at constant pressure is itself a function of temperature and can be approximated 
by an equation of the form 


Cpm = a+bT+cT^+dT 3 . 


(17) 


Molar entropy is Cpm/T so that 


Cpm/T = a/T +b+cT+dT 1 . 


(18) 
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Integrating with respect to T gives 

A Cpm= a ln(T/298)+b(T-29S)+ i fa c(7 2 -298 2 )+ 1 / 3 d(P-29 &) . (19) 

A surface describing AG for a chemical reaction, as a function of temperature and pressure, 
will show the regions of temperature and pressure in which the reaction is spontaneous. 


PROCEDURE 


The JANAF Thermochemical Tables 12 provide data for elements and compounds based on a 
reference state of 1 bar and 298.15 K. Values for heat capacity are given at a pressure of 1 bar in the 
temperature range of 0 to 6,000 K for most substances. Using an appropriate number of simultaneous 
equations, coefficients for best-fit curves were calculated from Cp values for species in the range of 
298.15 to 1,500 K. Species chosen were those included in equations (1) to (6). Best-fit curves were 
of the form 


Cpm = a+bT+cP 2 +dT ±3 . (20) 

Coefficients of T, as well as molar volume and standard Gibbs free energy change for species, are 
listed in table 2. 

A computer code was written to construct a three-dimensional surface of the change in Gibbs 
free energy of reactions (1) to (6) as functions of temperature and pressure. The generation of the 
surface is by equation (15) for solid species and equation (16) for gaseous species. Both equations 
(15) and (16) use equations of the form of equation (19) for the molar entropy term. In this program, 
temperature ranges from 298 to 1,500 K, and pressure ranges from 1 bar to 500 bars. The surface 
showing change in Gibbs free energy of a reaction subtracts LA G{T,p) of reactants from LAG(T,p) of 
products. 


RESULTS AND DISCUSSION 


Figures 1 a to 4 a show changes in Gibbs free energy of the reactions assumed relevant. The 
plots are surfaces as functions of temperature and pressure. Figure la is a plot of reaction 1: 

SiC+2H 2 = Si+CH4. (1) 

Except at lower temperatures and pressures, the reaction is spontaneous. Silicon is the stable solid, 
and CH4 is the gaseous product. The reaction is favored by increases in temperature and pressure. 
Initially, an increase in pressure decreases AG dramatically at all temperatures. At pressures 
greater than about 20 to 70 bars, the effect is less pronounced. Temperature has a significant effect 
at all pressures as well, and the effect is greater at higher pressures. This is consistent with Herbell 
et al. 1 Silicon carbide (SiC) should decompose at an increasing rate with increasing total pressure in 
that AG becomes more negative as pressure increases over the range of temperatures they 
considered. However, the rate is markedly less at higher pressures. Clearly, both temperature and 
pressure act in concert to enhance the decomposition of SiC in pure H 2 . 
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Figure 2 a shows reaction 2: 


SiC+H 2 0+H 2 = SiO+CHi. (2) 

A change in temperature has only a slight effect, and this reaction is not spontaneous at any 
temperature or pressure considered here. Silicon carbide is stable over the entire ranges of both 
parameters. It is interesting to speculate that the reaction is spontaneous at extremely high 
pressure, though this is extrapolation from the model. Should the trend continue, figure 2 suggests 
that, at extremely high pressures, SiC would convert to SiO and CH 4 in a hydrogen/water 
atmosphere. Solving for p in 

0 = AG°(SiO)+/?r ln(p)+AG°(CH 4 )+RT ln(p)-AG° (H 2 0)-/?7 ln(p ) 

- AG°(H 2 )-RT ln(p)-Vm(SiC)(p-l), 
the reaction would reach equilibrium at 40 kbars. 

The decomposition of Si 3 N 4 


Si 3 N 4 = 3Si+2N 2 (3) 

is shown in figure 3 a. This reaction is also not spontaneous at any temperature or pressure under 
consideration. In contrast to the previous reactions, this one is not favored by increasing temperature 
or pressure. In fact, figure 3 a suggests that the reaction will proceed only at very low temperatures 
and pressures. The pressure effect appears dominant at less than approximately 50 bars. Below 50 
bars AG, for the reaction, drops steeply and approaches zero only at much less than 1 bar. This is in 
keeping with Misra's 2 statement that reducing the pressure at constant temperature is equivalent to 
increasing the temperature at 1 atm. It is also a consequence of the fact that AG generally increases 
where the number of moles of product gases is greater than that of reactant gases. 

In the high-pressure range, temperature has a larger effect on AG than pressure. In this 
region, the change in entropy is of greater magnitude than the change in molar volume of 
constituents. 

Figure 4 a depicts reaction (4): 


Si 3 N 4 +6H 2 = 3SiH4+2N 2 . (4) 

The plot for this reaction is similar in appearance to those involving SiC, though products of this one 
are not stable over the ranges of temperature and pressure considered. Again, however, the reaction 
is favored by increasing temperature and pressure. 

Figure 5 a plots a possible reaction where Si0 2 forms on SiC: 

2SiC+2H 2 0+2H 2 = Si+Si0 2 +2CH4. (5) 
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The change in Gibbs free energy is the most negative of any reaction considered. It is therefore most 
likely to proceed. Except at low pressures, it is more dependent on temperature than pressure. As 
temperature increases, the change in entropy determines the reaction more than change in molar 
volume. Comparing figures 1 a and 5a, if SiC is unstable in pure H 2 , a wetter atmosphere provides 
some stability by allowing the formation of a protective Si0 2 scale. 

Figure 6 a shows a plot of the reaction forming Si0 2 on Si 3 N 4 : 

Si 3 N 4 +2H 2 0+3H 2 = Si+Si0 2 +SiH 4 +2NH 3 +N 2 . (6) 

Change in Gibbs free energy is not negative, so the reaction does not proceed as written. Instead, 
the opposite of the reaction proceeds. Again, Si 3 N 4 appears to be stable, though AG decreases with 
temperature. It is also apparent that AG for this reaction is the lowest for all reactions involving 
Si 3 N 4 in this study. This is most likely due to the great stability, hence negative AG of formation of 
Si0 2 over a broad range of temperatures. A positive AG for the reaction notwithstanding, it is less 
positive than for other reactions due to the possibility of Si0 2 as a product. 

Mathematical models do not generally reflect real processes completely accurately, yet 
models such as the one described here are illustrative and some tentative conclusions can be drawn. 
In a closed system, only SiC in pure H2 will decompose. Where moisture is present in the same 
proportion as hydrogen gas, SiC is stable up to moderate temperatures and pressures. Silicon nitride 
is stable in pure H 2 at all but the lowest pressures. This suggests that where temperatures are less 
than 1,500 K and pressures less than 500 bars, SiC will decompose quite readily in pure H 2 while 
Si 3 N 4 remains stable. Reliability of structural components and coatings in turbopumps and nozzles 
would be better served by using silicon nitride. 

A real structural material employs various sintering and binding aids which may engage in 
reactions with H 2 and H 2 0 where the parent material would not. Also, it can be argued that 
treatment of gases as ideal introduces error. Due to attraction between molecules, the molar volume 
of a real gas is less than that expected for an ideal gas. A smaller contribution by gaseous molar 
volume to overall Gibbs free energy might favor all of the above reactions as written. However, if the 
scale of energy change in reactions is large enough, the effect of additional terms may be 
immeasurably small. 

As an example, consider gases to behave as van der Waals gases, with molar volume 
derived in the following manner: 13 The compressibility factor is Z = pVm/RT for ideal gases. If 
p = RTI(Vm-b)-a/Vm 2 for a van der Waals gas, Z = Vm/(Vm-b)-a/RT. Dividing numerator and 
denominator by Vm gives 


Z = U{\-{b/Vm)]-a/VmRT . (21) 

At very low pressures, Vm is large and b/Vm is very small. A power series in Vm can be generated 
for the term \![\-{b/Vm)] to give 


Z= \+b/Vm+{b/Vm) 2 +...-dSVmRT 
= \+{b-a/RT){\IVm)+{b/Vm) 2 +... (22) 
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Since Z = pVm/RT, then Vm = ZRT/p = RT/p+RT/p(b-a/RT)(p/RTZ). More simply, 

Vm = RT/p+b-a/RT . (23) 


Integrating with respect to p gives 

Vm(p) = RT ln(p/pO)+(b-a/RT)(p-pO) . (24) 

If the units of the gas constant are /cJ/mol K and the units of (b-a/RT)(p-pQ) are m 3 bar/mol, the 
latter term must be multiplied by 100 to convert m 3 bar to kJ. 

Using equation (24) for molar volume in Gibbs free energy change for gases gives figures 1 b 
to 6b for the reactions investigated here. Table 3 gives van der Waals constants for gaseous 
species. 

In the range of temperature and pressure considered here, where AG is measured in kJ, there 
is no discemable difference in reactions where gases are treated as ideal or van der Waals gases. 
The implication is that gases behave ideally on the scale of energy changes considered here and can 
be modeled as such. At low pressure and/or high temperature, van der Waals attractive and 
repulsive forces have negligible effect, and a gas behaves very nearly ideal. However, where Gibbs 
free energy of reactions is calculated in kJ, the effect of van der Waals forces is simply too small to 
detect at temperatures and pressures investigated here. 


CONCLUSIONS 


On the basis of change in Gibbs free energy of reactions, Si 3 N 4 is more stable than SiC at all 
pressures and temperatures likely to be encountered in near-term advanced propulsion feed 
systems. Silicon carbide is more stable in a wet hydrogen atmosphere due to the formation of a Si02 
scale. Reactant and product gases can be treated as ideal for the purpose of reaction 
thermodynamics involving SiC, S^N,*, and H 2 /H 2 O. 
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Table 1. Reactions of SiC and Si 3 N 4 with H 2 and H 2 /H 2 O. 


High T. Low p 
Low T. High p 

High T. Low p 
Low T. High p 


Region I 

SiC+2H 2=Si+CH4 
SiC+2H2=Si+CH4 

Si3N4=3Si+2N2 

Si 3 N4 +6H 2 =3SiH4 +2N 2 


Regi on II 

S1C+H2 O+H 2=SiO+CH4 
SiC+2H2=Si+CH4 

Si3N4=3Si+2N2 
Si3N4+6H2=3SiH4 +2N2 


Region HI 


2SiC+2H20+2H2 =Si+SiO 2 
+2CH4 


Si 3 N4 +2H2 0+3H2 =Si02 + 
SH4 +2NH 3 +N 2 


1 , 

1 , 


Table 2. Thermochemical data. 


Cp Coefficients 


Species* 

AG g (kJ/moQ 

Vm 

a 

lz 

£ 

d 

h 2 (1) 

0.0 

- 

3.029xl0- 2 

-3.216xl0' 6 

3.102X10 9 

-2.045X10 4 

N 2 (1) 

0.0 

- 

2.487xl0- 2 

1.018xl0' 5 

-2.360xl0 9 

3.794xl0 4 

Si (3) 

0.0 

1.205xl0’ 5 

2.283x10-2 

3.839xl0‘ 6 

-3.536X10 2 

0.0 

SiC (2) 

-69.139 

1.249x1 O' 5 

-6.131xl0‘ 4 

1.174X10* 4 

-9.356x10'® 

2.578xl0- n 

S13N4 (1) 

-647.343 

4.412xl0‘ 5 

6.738x10-2 

1.241X10’ 4 

-3 .326x10-® 

-5.002X10 4 

SiO (2) 

-127.305 

- 

2.363x10-2 

2.520x10 - 5 

-1.746x10® 

4.357xl0- 12 

SiC>2 (3) 

-853.636 

2.590xl0' 5 

7.278x10-2 

1.285xl0‘ 6 

-4.141X10 3 

0.0 

SiH4 (1) 

56.827 

- 

2.107x10-2 

9. 134x10 - 5 

-2.781x10® 

-7.957X10 4 

H 2 0 (1) 

-228.582 

- 

2.823x10-2 

1.388xl0- 5 

-8.739xl0 10 

3.445xl0 4 

CH 4 (2) 

-50.768 

- 

2.229x10-2 

3.716xl0‘ 5 

3.095xl0' 3 

-1.81 lxlO" 1 1 

NH 3 (2) 

-16.367 

- 

2.634x10-2 

3.028X10’ 5 

4.560xl0' 9 

-4.581xl0- 12 


* Numbers in parentheses denote form of integrated equation used in entropy calculations. 

(1) : a ln(T) + bCT-TO) + dlin-TQ 2 ) - d/3(7'- 3 -7’0- 3 ) 

(2) : a ln(T) + b(T-T0) + c/KT^-TO 2 ) + miV-TW) 

(3) : a ln(T) + b(T-TO) + c/2(r 2 -r0' 2 ) 


Ref. 
2 , 8 , 11 
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Table 3. Van der Waals constants for gaseous species. 


Species 

aim 6 bar/mol 2 ! 

b (m 3 /mol) 

h 2 

2.476xl0" 7 

2.661xl0- 5 

n 2 

1.408x1 O ’ 6 

3.913xl0- 5 

SiO 

8.698xl0' 6 

8.582xl0' 6 

Sift* 

2.478xl0' 6 

3.275x10-5 

h 2 o 

5.536xl0- 6 

3.049x10-5 

CH 4 

2.283xl0' 6 

4.278x10-5 

nh 3 

4.225xl0' 6 

3.707x10-5 
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